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Rigid Body Mode Identification of the PAH-2 Helicopter 
Using the Eigensystem Realization Algorithm 

Axel Schenk and Richard S. Pappa 


Summary 

The rigid body modes of the PAH-2 'TIGER' helicopter were identified using the 
Eigensystem Realization Algorithm (ERA). This work complements ground vibration 
tests performed using DLR's traditional phase resonance technique and 1SSPA. Rigid 
body modal parameters arc important for ground resonance prediction. Time-domain 
data for ERA were obtained by inverse Fourier transformation of frequency response 
functions measured with stcppcd-sinc excitation. Mode purity (based on the Phase 
Resonance Criterion) was generally equal to or greater than corresponding results 
obtained in the ground vibration tests. All identified natural frequencies and mode 
shapes correlate well with corresponding ground vibration test results. 

The modal identification approach discussed in this report has become increasingly 
attractive in recent years due to the steadily declining cost and increased performance 
of scientific computers. As illustrated in this application modern time-domain meth- 
ods can be successfully applied to data acquired using DLR's existing test equipment. 
Some suggestions arc made for future applications of time domain modal identifica- 
tion in this manner. 


Acronyms 

APS Average Power Spectrum (equation 1) 

ATWL Analysis Time Window Length 

CMI Consistent Mode Indicator (ERA'S primary accuracy indicator) 

DLR Deutsche Forschungsanstalt fur Luft- und Raumfahrt e.V. 

EMAC Extended Modal Amplitude Coherence 

ERA Eigensystem Realization Algorithm 

FRF Frequency Response Function 

GVT Ground Vibration Test 

IRF Impulse Response Function 

ISSPA Identification of Structural System Parameters (method) 

MAC Modal Assurance Criterion (mode shape correlation measure) 

(equation 4) 

MBB Messerschmitt-Bolkow-Blohm GmbH, a division of Deutsche Aero- 

space 

MIF Mode Indicator Function (equation 2) 

PAH Panzer Abwehr Hubschrauber (anti-tank helicopter) 

PRC Phase Resonance Criterion (equation 3) 

PRM Phase Resonance Method 




1. Introduction 


In January 1991, ground vibration tests were performed by DLR on the first proto- 
type of the Franco-German PAH-2 'TIGER' helicopter [1] at the MBB site in Otto- 
brunn, Germany [2]. Test objectives were identification of structural dynamic 
characteristics (modal parameters) and measurement of transmissibilities (fre- 
quency response functions). Among several issues of helicopter dynamics, the 
ground resonance problem is of great significance. Ground resonance is a 
dynamic instability involving the coupling of the rotor blade lag motion with the 
in-plane motion of the rotor hub [3], When in contact with the ground, the critical 
structural mode is usually a rigid body oscillation of the helicopter on its landing 
gear. Ground resonance instability is avoided by having sufficiently-low rigid 
body natural frequencies and/or adequate damping. 

As a complement to the ground vibration tests conducted using DLR's traditional 
phase resonance technique [4] and ISSPA [5,6], some post-test time domain 
analyses were also performed in August 1991. The rigid body modes were 
selected as the target of this work because frequency response functions were 
available from [2], and as a contribution to the prediction of ground resonance 
behavior. Time domain data (specifically, impulse response functions) were 
derived from the measured frequency response functions using an inverse fast 
Fourier transformation. 

In this report the application and results of the time domain analysis are pre- 
sented. The Eigensystem Realization Algorithm (ERA) [7-9] was used for modal 
parameter identification. Recently, a comparison of results obtained with the 
phase resonance method (sinusoidal excitation with appropriated forces) and with 
ERA (multiple-input analysis of directly measured free-decay time histories) has 
shown that these dissimilar methods both yield reliable modal parameters of high 
correlation [10]. Correlation of results provides increased confidence in difficult 
testing situations. 

The modal identification approach discussed in this report has become increas- 
ingly attractive in recent years due to the steadily declining cost and increased 
performance of scientific computers. As illustrated in this application, modern 
time-domain methods can be successfully applied to data acquired using DLR's 
existing large test facility. No additional tests are necessary by using frequency 
response functions acquired during initial sweep runs. 


2. Test Article 


Figure 1 shows a photograph of the test setup. During the tests the helicopter was 
standing on the ground and was configured as follows: 

• no main rotor blades, 

• no tail rotor blades, 

• no mast mounted sight, 
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• with dummy masses for horizontal stabilizer end plates, 

• with normal shock absorber pressure, 

• with normal tire pressure, 

• with pilot and gunner dummy masses (2 x 110 kg), 

• with fuel (1088 kg), and 

• with dummy stubwing stores (480 kg). 


The gross weight during testing was approximately 4600 kg and was measured 
means of a force transducer. The overall dimensions of the structure are: 


by 


• length : 15.03 m, 


• height : 5.09 m, 

• width : 4.34 m. 


3. Data Acquisition and Assessment 


F.ve sets of frequency response functions (FRFs) were available for analysis The 

fnZ?nn!Z k 6a ^ t Set are listed in Tab,e The corresponding exciter 
locat.ons are shown in Figure 2. These FRF data were measured using stepped- 

n^nn an , a i0n c ° nve 4 rted ,nto di 9 ital| y filtered (bandlimited) time histories 
funrUnn^ hrfT T* t , ra ^ forma,ion - Each set of resulting impulse response 
frrm.pnrx RFS CO " S ,i S S of 242 acceleration histories. Because of the inconsistent 
frequency ranges and number of data points of the various tests, each data set is 
analyzed separately in a single-input, multiple-output analysis. 

Representative FRFs and corresponding IRFs are shown In Figures 3 through 7 
These data were measured on the horizontal stabilizer in the x, y, and z directions' 
respectively, going from top to bottom in each figure. The noise level is very low 
based on the smoothness of the FRF curves. Low noise, of course, simplifies data 
analysis. However, several characteristics which complicate data analysis are 
also seen. First, nonlinear behavior (unsymmetrical response peaks) is observed 

near 2 4 Hz in Record 183 of Figure 4. This observation 
correlates with the discussion of nonlinearities in [2], associated primarily with the 
achons of the wheels and landing gear. Secondly, spurious data spikes occur near 
4.9 Hz in all three FRFs in Figure 6. These "gain peaks" occurred when sensor 
amplifiers were switched to avoid overloads near resonance. Thirdly more than 
six response peaks occur in some FRFs, for example, in Record 183 of Fiqure 3 
However, only the six rigid body modes are known to occur in this frequency 
interval. A fourth data irregularity is observed in Record 183 of Figure 3 in the 
frequency interval from 4.6 Hz to 6.0 Hz. No explanation can be given for the 
unusually-flat, large magnitude function seen here. 


The transformation of bandlimited FRFs into IRFs calculates a fictitious sampling 
frequency based on the frequency range of the FRFs. Therefore, the data oscil- 
lation frequencies observed in the IRFs of Figures 3 through 7 are distorted The 
damping rates, however, are unaffected by the transformation. Using Record 181 
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of Figure 3 as an example, this distortion is explained as follows. The first peak in 
the FRF occurs at approximately 0.25 Hz after the initially measured frequency of 
1.4 Hz. The associated cycle time is computed as the reciprocal of 0.25 Hz or 
4.0 sec. This cycle time is observed in the IRF for the dominant oscillation. How- 
ever, the correct frequency for this peak is, in fact, 1.65 Hz, corresponding to a 
cycle time of 0.61 sec. During analysis, this distorsion is eliminated by adding the 
initially measured frequency to the frequency results. This is a standard feature 
of ERA. With bandlimited FRFs, the sampling frequency of the IRF equals twice the 
bandwidth of the corresponding FRF. Similarly, the number of time samples in the 
IRF equals twice the number of frequency lines in the FRF. 

As an overview of all 242 FRF measurements, Figures 8 through 12 show average 
power spectra (APS) and mode indicator functions (MIFs) for each of the five data 
sets. These functions are calculated as follows: 


( 1 ) 

and 


APS (0 = 


El M0 I 2 


( 2 ) 


MIF (f) = 


1.0 


E I"/ (01 I M0 I 

El MO I 2 


• 1000 


where u t (f) and | u, (f) | are the real part and magnitude, respectively, of the /th 
response measurement with respect to the applied force. n p is the total number of 
physical response measurements. 


Both of these functions display peaks at each natural frequency. Additionally, the 
APS shows the relative magnitude of each modal response. The MIF, on the other 
hand, provides no information concerning modal response magnitudes; however, 
the resolution of individual modes is much higher. As noted previously in the FRF 
plots, significant deviations from linear behavior (unsymmetrical response peaks) 
are seen in some functions, for example, for the mode near 2.5 Hz in the MIF of 
Figure 9. Also, the "gain peak" observed in test AFXZL1 is clearly seen as a large 
spike in the APS of Figure 11. 

In summary, appreciable nonlinearity and other minor data irregularities are 
observed in the FRFs measured on the PAH-2 helicopter. Such characteristics are 
not uncommon in dynamics tests of large, built-up structures. These factors com- 
plicate, but do not preclude, accurate modal identification. 
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4. Initial Results 


4-1 ERA Application Approach and Evaluation Criteria 

Previous experience has shown that significant changes can occur in ERA ana- 
lyses as a function of the assumed number of modes. In particular optimum 
accuracy for different modes typically occurs at different numbers of assumed 
modes Also weakly excited modes often require relatively high numbers of 
assumed modes to be properly identified. For these reasons, the assumed number 
of modes is incremented in each analysis from 1 up to 15. 

Reliab'e modal parameters, in particular the damping values, are obtained if 
ERA s primary accuracy indicator, the Consistent-Mode Indicator (CMI) reaches 
values larger than approximately 80 percent. Compared with the damping results 
identif'ed frequencies and mode shapes are typically reliable at lower CMI values’ 
CMI measures the consistency of identified mode shape components both tem- 

a re ^ o ni Da red a \l/M 1!^ C ° mp0ne H nts from the beginning of the da?a analysis window 
are compared with corresponding components at the end of the window Addi- 

calculation 6 CMMs col,ln ear.ty of mode-shape phase angles is included in the 
is givenTn [11] Pr ° Ven ’ accuracy indicator. A complete description 

terh^lpRC^ 6 PRC^ is' nnt S a eS f Se< H a 6 wel,establls hed Phase Resonance Cri- 
erion (PRC). PRC is not a standard ERA accuracy indicator, but is used here for 

compansons wilh DLR s tradilional phase resonance testing approach I. ifcom- 
puled in the same manner as MIF, Eq. (2), except using the identified (complex) 
mode shapes instead of response measurements: P 


( 3 ) 


PRCj = 


H 

£ i % 1 1 % 


1.0 


/= i 


Ei*. 


/= i 


1000 


where <D,' and | <D„ | are the real part and magnitude, respectively of the Ah 
component of the;th identified mode shape. Modes with PRC values greater than 
950 are generally accepted as well-identified normal modes. 


4.2 Frequency, Damping, Accuracy Indicators, and MAC 

Initial ERA results for each of the five tests are shown in Figures 13a throuah 
17a, respectively. Each figure consists of a frequency plot, a plot of the viscous 
damping factor (percent of critical damping), a plot of CMI (as a measure of reli- 
abi ity) and a plot of PRC (as a measure of normal mode purity). These results are 
plotted versus the assumed number of modes, varying from 1 to 15. Each row of 
results corresponds to a separate ERA analysis with a specified number of 
assumed modes. 
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In the frequency plot, each detected mode is represented by a vertical dash at the 
associated frequency. The confidence in each result is expressed by the length of 
the vertical dash which is proportional to CMI. The highest confidence is attained 
if the distance between minor tic marks on the vertical axis is filled. Thus, more 
confidence is associated with results appearing as solid vertical lines than with 
those appearing as dotted or dashed lines. 

In the remaining three plots in each figure, the corresponding results for damping, 
CMI, and PRC are displayed using numerals for each mode appearing in the fre- 
quency plot. The parameters for the lowest frequency mode (along each row) are 
plotted using a "1", for the second-lowest frequency (along the same row) using a 
"2", and so forth. Due to the actual number of identified frequencies, the number- 
ing of damping values and accuracy indicators of individual modes may change 
over successive rows. 

The identified frequencies in Figures 13a through 17a are labeled with the defor- 
mation pattern based on associated mode shape results. Using the Modal Assur- 
ance Criterion (MAC) [12], each of the ERA-identified mode shapes is correlated 
with each of the six rigid body modes determined in the ground vibration test 
[2], as follows: 

I I 

( ^ GVT ^GVT ) ( ^ERA ®ERA ) 

where H designates the Hermitian transpose. MAC values range from 0.0 for two 
independent (orthogonal) vectors to 1.0 for two identical vectors (differing at most 
by a scale factor). Values greater than approximately 0.7 indicate a considerable 
degree of similarity. 

These MAC results are shown in Figures 13b through 17b using a simple graphical 
format. Each row and column in these figures represents one mode shape. The 
value of MAC is proportional to the size of the rectangle drawn at the intersection 
of the corresponding row and column. MAC values greater than 0.7 are darkened 
for emphasis. In general, each ERA-identified mode correlates highly and uniquely 
with one of the six modes identified in the GVT. Each mode is labeled with three 
numbers. The first of which is the case number, i.e. the number of assumed 
modes. The second of which is the order number of the mode in the particular 
case. The associated frequency is also shown as the third number. Best mode 
shape results based on PRC from the GVT are recorded along the ordinate while 
all ERA mode shapes are recorded along the abscissa. 

All GVT modes are decoupled from each other. Thus, mode coupling of ERA 
results is indicated if in the same column one ERA mode correlates with two or 
more GVT modes, for example in Figure 16b, case 15: ERA mode 4 with GVT 
modes 4 and 6. On the other hand, repeated identification of the same mode shape 
within the same case is also indicated if several ERA modes correlate with one 
GVT mode, for example in Figure 15b, case 15: ERA modes 1, 3, and 4 with GVT 
mode 1. Ideally, each rigid body mode should be identified only once in each case, 
for example in Figure 15b, yaw mode. 


( 4 ) MAC GVT I ERA ~ 
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The initial results obtained for each test are discussed individually in the following 
paragraphs: 

For test AFXL1, Figure 13, stable results in both frequency and damping are 
achieved for yaw, vertical, and roll modes. Stable frequency results are also 
obtained for the longitudinal mode; however, the corresponding damping factors 
show greater spread at values between 6 and 8 percent. Increasing the number 
of assumed modes, the pitch mode damping decreases continuously to values 
around 1 percent while CMI increases to approximately 50 percent. The pitch 
mode also shows "mode splitting" at three particular values of assumed number 
of modes (10, 14, and 15), attributed to nonlinearity. Note in Figure 13b that the 
corresponding pairs of mode shapes all show high correlation with the GVT pitch 
mode. With the exception of the pitch mode, best CMI results for all detected 
modes are higher than 80 percent and PRC values are higher than 950. 

Test AFXL2, Figure 14, used the same exciter location as test AFXL1 (longitudinal 
excitation at the rotor hub) except at a higher force level (160 N vs. 140 N). Again, 
stable frequency and damping results are obtained for yaw, vertical, and roll 
modes. Due to nonlinearity, the frequencies are 25 percent (yaw), 12 percent (ver- 
tical), and 7 percent (roll) higher than the corresponding frequencies of test AFXL1. 
Now, best CMI values of these modes have increased to above 90 percent and 
PRC values to above 970. Based on mode shape correlations (Figure 14b), multi- 
ple identification of the longitudinal mode occurred in the frequency range of 
approximately 2 to 3 Hz. The associated damping estimates show wide variation 
between 2 and 7 percent. Best CMI values are approximately 80 percent. Also, 
PRC values higher than 950 (associated with the 2.35 Hz results) indicate good 
mode purity. Considering the highly unsymmetrical peak in the MIF and the high, 
but flat-top response in the average power spectum (Figure 9), the unusual phe- 
nomenon of multiple identification may be explained and attributed to nonlinearity. 
Since the measurements were limited to frequencies below 5 Hz, no pitch mode is 
identified in test AFXL2. The lateral mode is not detected in test AFXL2 (or AFXL1) 
because of the inappropriate excitation direction. 

Although test AFYL1 was performed to excite the lateral mode better (using lateral 
excitation at the rotor hub), surprisingly only one acceptable identification of this 
mode is obtained (at 14 assumed modes), see Figure 15. However, a CMI value 
of only 48 percent indicates relatively low confidence. At 15 assumed modes, 
multiple identification occurred with even lower CMI than in the previous case! 
Since this mode is identified all together only five times (in three cases), an eval- 
uation of stability of modal parameters is impossible. The CMI results for the lon- 
gitudinal and vertical modes are similarly low due to the inappropriate direction 
of excitation. On the other hand, stable frequency and damping results are 
obtained for both the yaw and roll modes. These modes attain CMI values above 
80 percent and PRC values above 970. 

The mode indicator function of test AFXZL1, Figure 11, shows several peaks below 
3 Hz. Based on mode shape correlations, modes are extracted in this frequency 
range which are combinations of lateral and longitudinal, and vertical motions, 
coupled with pitch, see Figure 16. The corresponding damping results show a wide 
spread of values between 6 and 10 percent. Also, CMI as well as PRC indicate low 
confidence in these modes. Once again, yaw and roll modes are identified with 
stable frequency and damping values, and CMI values of approximately 80 percent 
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and PRC values above 950. The pitch mode is identified at 4.9 Hz with a best CMI 
value of 56 percent and a PRC value of 929. Note that this frequency coincides 
exactly with the "gain peak" discussed in Section 3. However, the plot of the mode 
indicator function, Figure 11, discloses that a response peak is at exactly 5 Hz 
which can be attributed to the pitch mode. Since the gain peak generally extends 
higher in magnitude than the true FRFs, an incorrect frequency result is apparently 
obtained for the pitch mode. In addition, the damping result decreases to zero with 
increasing assumed number of modes which is in good agreement with a narrow 
peak in the FRFs. Obviously, both frequency and damping results for this mode 
are questionable. On the other hand, high confidence can be given in the mode 
shape estimate because a pitch mode pattern is obtained. 

Figure 17 shows the initial ERA results for test AFXZL2, using antisymmetric exci- 
tation at the landing gear. Low-confidence results with CMI less than 40 percent 
are again found for both lateral and longitudinal modes. CMI values increase to 
95 percent for the yaw mode with increasing numbers of assumed modes, while for 
the roll mode CMI values are optimum at moderate numbers of assumed modes. 
Nevertheless, stable frequency and damping estimates and PRC values above 960 
are obtained for both modes. A mode of low confidence (with CMI of approxi- 
mately 20 percent and PRC near 800) is found at a frequency of 5.1 Hz which cor- 
relates with none of the six rigid body modes. Although the CMI values are not 
quite satisfactory, the pitch mode is clearly detected at 5.6 Hz, with damping and 
mode shape based on PRC showing stable results. 

The best results for each mode based on CMI are listed in Table 2, together with 
the analysis time window length (ATWL) used in each case. Only certain modes 
are identified in certain tests. Also, most modes show a significant spread of fre- 
quency and damping values. This can be explained by the different excitation 
locations and levels used in the various tests and the nonlinear response of the 
structure attributed to the actions of wheels and landing gear. 

In summary, four of the six rigid body modes are identified with high confidence 
(CMI of at least 80 percent) in the initial ERA analyses. The remaining two, the 
lateral and pitch modes, can only be obtained with unsatisfactory confidence. 
Further analyses are needed, in particular to evaluate reliable damping factors. 
Appreciable variations of modal parameters (frequency and damping) occurred for 
all six modes due to structural nonlinearities. Comparing the FRF measurements, 
Figures 3 through 7, with the ERA findings, Figures 13 through 17, good agreement 
in corresponding frequencies is seen. 


4.3 Mode Shapes 

Identified mode shapes from the initial analyses are plotted in Figures 18 through 
23. Modes having the highest CMI values in Table 2 are shown. All six modes 
correlate to a high degree with the corresponding mode shapes identified in the 
GVT using the traditional phase resonance method and ISSPA. Correlation coeffi- 
cients (MAC) of the GVT and initial ERA results are listed in Table 3. 
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5. Improved Results 


The motivation for the improvement of results, discussed in this section, can be 
seen in Table 2. Most of the modes have CMI values of approximately 8o’ percent 
or higher and sufficient confidence can be placed in the associated modal param- 
eters. On the other hand, some modes show insufficient CMI values. This low 
confidence in some results is also expressed by dotted or dashed frequency lines 
in Figures 13a through 17a, for example in Figures 15a through 17a for the lateral 
mode. Together with the longitudinal mode, the lateral mode is most important to 
the ground resonance problem. Therefore, accurate frequency and damping esti- 
mates for these modes are a major objective of the test program and, also, of the 
ERA analyses reported here. 


5.1 Approach 

Complex structures require a high number of sensors for adequate resolution and 
identification of all mode shapes. However, only a subset of all sensors attached 
to the structure is needed for identification of specific modes. Modal identification 
of target modes, in particular the damping estimation, can be improved if a subset 
of the inputs (excitation locations) and outputs (sensor locations) is emphasized 
which matches the requirements of the target modes. Then, unsuitable excitations 
and sensors with lower signal-to-noise ratios are neglected. While emphasizing 
the target mode, normally all other modes are de-emphasized. This secondary 
feature fails if the selection of Key excitations and/or sensors accidently matches 
other modes. Each mode of interest normally requires a dedicated analysis. Suc- 
cessful applications of this improvement technique have been demonstrated in 
[9.10]. 

In this project, mode shape patterns for selection of key sensors are provided by 
the initial analyses. Ten optimal responses are used in each improvement analy- 
sis. Criteria for selection are high accuracy indicator values (EMAC) and a corre- 
sponding high modal amplitude at the location of the key sensor. EMAC, the 
Extended Modal Amplitude Coherence, is used in the calculation of CMI and pro- 
vides information concerning the controllability and observability of the modes in 
terms of each individual excitation point (input) and sensor measurement (output). 

A feature of ERA is that complete mode shapes are derived in analyses using key 
sensors. This is accomplished by retaining all response measurements in the first 
block row of ERA'S data input matrix. In all following block rows, only the 
responses of the specified key sensors are included. 


5.2 Results 

Six results from Table 2 are selected as target modes for the improvement tech- 
nique described above. CMI values of 

- lateral mode Tests AFYL1 and AFXZL2; 

- vertical mode Test AFYL1; 
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- pitch mode Tests AFXL1 , AFXZL1 , and AFXZL2 

are below 70 percent. The associated damping values show some variation and, 
in particular for the lateral mode, only a few estimates were obtained in the initial 
analyses. 


5.2.1 Lateral Mode 

Improved results for the lateral mode are shown in Figures 24 and 25 as expanded 
views of the frequency range from 1.0 to 2.0 Hz. Best CMI values in both analyses 
are now approximately 80 percent and PRC values in test AFYL1 have increased 
considerably to above 970. In both damping plots, some random variation is seen. 
However, a damping value between 6 and 12 percent is estimated with the excita- 
tion used in test AFYL1 (also in other analyses not shown here), while a damping 
between 10 and 15 percent is estimated with the excitation used in test AFXZL2. 

The response amplitude of the lateral mode was approximately twice as large in 
test AFXZL2 than in test AFYL1 , based on [2]. Due to nonlinearities, the natural 
frequency of this mode decreased from 1.44 Hz in test AFYL1 to 1.23 Hz in test 
AFXZL2, with a corresponding damping increase by a factor of two. These nonlin- 
earity trends correlate closely with those reported in [2], 

5.2.2 Vertical Mode 

The current limitation of the improvement technique is demonstrated with the ver- 
tical mode in test AFYL1, Figure 26. Little improvement is possible due to the 
inadequacy of the excitation in this test. The best CMI value obtained is now only 
56 percent and the associated PRC value is 883. This is an example where two 
other strongly excited modes in the vincinity of the target mode are not signif- 
icantly deemphasized. The lateral excitation used in this test is more appropriate 
to the yaw and roll modes than to the vertical mode. 


5.2.3 Pitch Mode 

Expanded views of the frequency range from 4.5 to 5.8 I- z are provided in Figures 
27 through 29, showing improved results for the pitch node. The wide variation 
of frequency and damping according to different excital ons used in tests AFXL1 , 
AFXZL1 and AFXZL2 and seen in Table 2 is corroborated by the improvement 
results. 

A good stabilization of damping results with high CMI values is achieved in tests 
AFXL1 and AFXZL2, while no stabilization is obtained in test AFXZL1. In the latter 
case, results with CMI above 50 percent have damp ng values between 1.0 and 
2.6 percent. The low CMI values of AFXL1/14 and of AFXZL1/12 in the initial ana- 
lyses are increased to values of 80 percent or higher together with some increase 
of PRC. Smaller accuracy improvement from 68 to 77 percent in CMI is possible in 
test AFXZL2; however, the corresponding PRC valu 3 decreases from 966 to 941. 
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5.3 Assessment 


S.Hinn« « ° f Tab 6 2 We ? , sel ® cted as candidates for improvement analyses. The 
findings are summarized in Tables 4 and 5. With the exception of the vertical 
mode, significant improvement of accuracy based on the parameter CMI was 
achieved Compared to the results shown in Table 2, more confidence can now be 
given in the stated modal parameters. In five of the six cases, CMI values have 
increased considerably. With the vertical mode, however, less improvement of 
CMI was achievable due to an insufficient excitation and other dominant modes 
at nearby frequencies. 


6. Comparison of ERA and GVT Results 


A comparison of ERA and GVT results is given in Table 6. Only modes with similar 
exci ation are compared. Unfortunately, no damping estimates were available from 
12J for these modes. Therefore, best identified modes based on PRC from the 
ground vibration tests are also listed and indicated with an asterisk. 

Overall, there is excellent correlation of frequency results. As discussed earlier 

I/«h,^ 0rr f SP °K dln9 m 5 de shapes are also h '9 h, y similar. The identified damping 

fh l ® J 1 ? 0 Sh ° W f ?° od a 9 reemen t although, of course, the sparsity of results limits 
the usefulness of the comparison. 

PRC values obtained using ERA are higher in 9 of the 10 cases than correspondino 

the U GVT° b amed f h th Z ? VT ' , M mUS< be realized that a tj 9ht time schedule P during 
the GVT prevented additional efforts to achieve higher PRC values 9 


7. Conclusions 


Rigid body mode identification of the PAH-2 'TIGER' helicopter was performed 

orm,nri h v e ih E ^ enS r l ? m Raal ' zat ' on Al 3 0rilhm (ERA). This work complements the 
S™", *' bza,on '“Is performed with DLR s classical PRM approach and ISSPA. 

re^abll wSh pbm Z SilT I 0 ' ! he dampin9 values ' could ba svaluated 
renaoiy with PRM and ISSPA due to nonlinearities. 

fpa f fS , WHh various J orce level * and excitation locations were analyzed 
2 ^RA. Analyses were performed in a single-input manner due to differences in 
data acquisition parameters among the data sets. Nevertheless, reliable results for 
most modes were obtained in the initial analyses. A wide variation of modal 
parameters was found among the separate tests due to the highly nonlinear 
response of the structure. For modes with low accuracy indicators 9 improvement 

fota! y numb^ r of S 242 : . eSSjU,,y C ° ndUC,ed em P hasizi ng a subset of ten sensors of a 

All identified frequencies and mode shapes correlate well with their countemarts 

n »»L 9 r nd Vibrali0a ^ 1651 [2] -, PRC va,uas ° blai "«l With ERA we e equaMo or 
greater than corresponding results obtained in the GVT. 
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In closing, the authors would like to make the following suggestions for improving 

future applications of ERA: 

• Perform ERA analyses at the test site using FRF data measured during initial 
PRM sweep runs. These sweep runs are always conducted prior to tuning 
individual modes; however, FRFs are not normally saved during this phase 
of testing. Only the mode indicator function is usually evaluated. ERA-identi- 
fied mode shapes would be available to assist in subsequent selection of 
optimum shaker locations and force distributions. 

• Use the same setup of data acquisition parameters for all tests in order to 
permit multiple-input ERA analyses. In this case, further reduction of analysis 
time is possible. With sequential single-input analyses, mode shape corre- 
lation is needed to associate results for the same mode. This step is unnec- 
essary with multiple-input analysis because only one result per mode is 
extracted as the best least-squares estimate. With highly nonlinear responses, 
however, multiple-input analysis may reach its limitation. If input forces and 
exciter locations used in the individual tests cause too great a variation of 
structural response and accordingly shifts in modal parameters, single-input 
analysis may be preferable (or multiple-input analysis with a reduced number 
of selected inputs). 

• Begin sweep runs at a frequency sufficiently less than the lowest frequency 
of interest. Sweep runs which begin too close to a mode of interest, for 
example, the mode at 1.5 Hz in Record 182 of Figure 3, cause too great a 
distortion on the corresponding IRF for accurate modal identification, in par- 
ticular for the damping value. 

• Minimize data acquisition time by using unequally spaced frequency lines. A 
minimum of five data points should be measured in the half-power bandwidth 
of each resonance, but the number of data points acquired between reso- 
nances can be much less. 
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Test 

Excitation 
Loc./Force 
- / N 

Frequency 

Range 

Hz 

Time 

Samples 

Time 

Window 

sec 

Sampling 

Frequency 

Hz 

AFXL1 

AFXL2 

AFYL1 

AFXZL1 

AFXZL2 

El / 140 
El / 160 
E2 / 120 
E21 ,E22 / 500 
E21.E22 / ±500 

1.4 - 6.0 
1.3 - 5.0 
1.0 - 5.0 
1.0 - 6.0 
1.0 - 6.0 

184 

370 

400 

500 

200 

20.0 

50.0 

50.0 

50.0 

20.0 

9.2 

7.4 

8.0 

10.0 

10.0 


Table 1. Test Parameters (Sine Sweeps) 


Mode 

Test/Case 

Frequency 

Damping Factoi 

CMI 

PRC 

ATWL 



Hz 

(viscous), % 

% 


sec 

Lateral 

AFYL1/14 

AFXZL2/10 

1.22 

1.31 

4.50 

11.71 

48 

36 

934 

963 

4.457 

4.100 

Longitudinal 

AFXL1/09 

AFXL2/15 

1.65 

2.35 

6.60 

2.99 

90 

78 

960 

976 

4.457 

5.541 

Yaw 

AFXL1/14 

AFXL2/14 

AFYL1/12 

AFXZL1/15 

AFXZL2/14 

2.78 

3.48 

3.41 

3.53 

3.46 

1.78 

2.27 

2.16 

2.42 

2.43 

85 

95 

88 

92 

95 

966 

973 

972 

960 

968 

4.457 

5.541 

5.125 

4.100 

4.100 

Vertical 

AFXL1/1 1 
AFXL2/14 
AFYL1/13 

3.38 

3.77 

3.81 

2.39 

2.26 

2.30 

83 

93 

53 

951 

986 

861 

4.457 

5.541 

5.125 

Roll 

AFXL1/15 
AFXL2/1 1 
AFYL1/13 
AFXZL1/12 
AFXZL2/09 

4.00 

4.26 

3.96 

4.02 

3.82 

1.11 

1.06 

1.86 

3.12 

4.94 

94 

96 

89 

79 

83 

977 

979 

971 

972 
968 

4.457 

5.541 

5.125 

4.100 

4.100 

Pitch 

AFXL1/14 

AFXZL1/12 

AFXZL2/13 

4.59 
4.90 

5.59 

0.91 

1.72 

2.57 

54 

56 

68 

928 

929 
966 

4.457 

4.100 

4.100 


on CMI) 


Table 2. Initial Identification Results (Best Results Based 




Mode 

Mode Reference Label 

MAC, 


GVT 

ERA 

% 

Lateral 

AY101 

AFYL1/14 

99 

Longitudinal 

A09L2 

AFXL1/09 

91 

Yaw 

A05L1 

AFXL2/14 

97 

Vertical 

A03L1 

AFXL2/14 

90 

Roll 

AY103 

AFXL2/1 1 

96 

Pitch 

A02L1 

AFXZL2/13 

87 


Table 3. Correlation of GVT and Initial ERA Mode Shapes 


Mode 

Test/Case 

Frequency 

Hz 

Damping Factor 
(viscous), % 

CMI 

% 

PRC 

ATWL 

sec 

Lateral 

AFYL1-KD4/10 

1.44 

7.79 

78 

974 

2.875 


AFXZL2-KD3/14 

1.23 

12.55 

78 

964 

4.100 

Vertical 

AFYL1-KD1/14 

3.83 

2.94 

56 

883 

5.250 

Pitch 

AFXL1-KD5/1 1 

4.61 

0.96 

85 

959 

7.935 


AFXZL1-KD1/12 

4.90 

1.96 

80 

943 

4.100 


AFXZL2-KD4/13 

5.61 

2.60 

77 

941 

4.100 


Table 4. Improved Identification Results 


Mode 

Mode Reference Label 

MAC, 


GVT 

ERA 

% 

Lateral 

AY101 

AFYL1-KD4/10 

99 

Pitch 

A02L1 

AFXL1-KD5/1 1 

93 


Table 5. Correlation of GVT and Improved ERA Mode Shaf es 
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1C 

1 ERA 

974 

964 

960 

976 

972 

968 

986 

968 

971 

943 


PF 

GVT 

976 

930 

940 

930 

958 

803 

657 

971 



840 

927 

Damping 

:tor, 

ERA 

7.79 

12.55 

6.60 

2.99 

2.16 

2.43 

2.26 

4.94 

1.86 

1.96 


Viscous 

Fac 

0 

GVT 

1 1 

3.38 $ 

2.18 $ 

1.83 $ 

l 1 

2.43 $ 


lency, 

Iz 

ERA 

1.44 

1.23 

1.65 

2.35 

3.41 

3.46 

3.77 

3.82 

3.96 

4.90 

indicates best GVT result based on PRC in [2] 

$ indicates average value of several damping estimates 

Freqt 

K 

GVT 

1.46 

1.30 

1.64 

2.42 

1.77 

3.44 

3.45 
2.52 

3.78 

3.80 

O CO 
CO CTj 

CO CO 

4.98 

5.22 

Reference Label 
ERA 

AFYL1-KD4/10 

AFXZL2-KD3/14 

AFXL1/09 
AFXL2/15 
no FRFs available 

AFYL1/12 
AFXZL2/14 
no FRFs available 

AFXL2/14 
no FRFs available 

AFXZL2/09 

AFYL1/13 

AFXZL1-KD1/12 
no FRFs available 

Mode 

GVT 

AY101 * 
AXZ201 

A10L1 
AX201 
A09L2 * 

AY102 
AXZ202 
A05L1 * 

AX202 
A03L1 * 

AX2203 
AY103 ’ 

r\J * 
o 

N cd 

>> o 
< < 

Mode 

Lateral 

Longitudinal 

Yaw 

Vertical 

Roll 

\5 

I 


Table 6. Comparison of ERA Results with GVT Results from [2] 


indicates best GVT result based on PRC in [2] 
indicates average value of several damping estimates 
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Figure 2. Exciter Locations 
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Figure 13a. Initial ERA Results for Test AFXL1 
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Figure 13b. Correlation of ERA and GVT Modes for Test AFXL1 
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Figure 14b. Correlation of ERA and GVT Modes for Test AFXL2 
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Figure 15a. Initial ERA Results for Test AFYL1 










MODAL ASSURANCE CRITERION (MAC) 


- 33 - 



*96*£ 

II 

/Si 

60B‘£ 

01 

/si 

eike 

8 

/SI 

916*2 

9 

/SI 

0U*2 

t 

/Si 

l£2*l 

E 

/Si 

2£2*l 

l 

/si 

E96E 

Ot 

Ai 

£18*£ 

B 

/m 

£ IKE 

L 

/M 

££6*2 

S 

/n 

tZZ'i 

l 

Ai 

E96*E 

6 

/si 

2 18*£ 

B 

/El 

lift 

9 

/ti 

9S62 

S 

/Cl 

l£2* l 

l 

/Cl 

t96*E 

8 

fz l 

S18*£ 

L 

/2l 

2 IKE 

6 

/2l 

S962 

t 

/2l 

6tS> 

6 

/u 

996‘E 

2 

/u 

2 18*£ 

S 

/n 

9 IKE 

► 

/u 

S26*2 

E 

/u 

66*> 

B 

/Ol 

£26*£ 

9 

/Ol 

OIB’E 

6 

/Ol 

BIKE 

t 

/Ol 

£26*2 

E 

/Ol 


B26'£ 9 /6 

C2B*£ S /6 

8 IKE t /6 

226*£ 9 /8 

lEB'E S /B 

2 IKE » /B 

226*£ 9 /L 

tU£ S /L 

12KE £ /Z 

126*E * /9 

t08*£ £ /9 

»2K£ 3 /9 

286*E E /« 

22KE 2 /£ 

9CKE Z /£ 


(iao woaj s3Q0h) ANrvs3aow~nad 


Figure 15b. Correlation of ERA and GVT Modes for Test AFYL1 
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Figure 16a. initial ERA Results for Test AFXZL1 
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Figure 16b. Correlation of ERA and GVT Modes for Test AFXZL1 
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Figure 17b. Correlation of ERA and GVT Modes for Test AFXZL2 
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Figure 18. Lateral Mode Shape 
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Figure 19. Longitudinal Mode Shape 
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Figure 21. Vertical Mode Shape 
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Figure 23. Pitch Mode Shape 
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Figure 28. improved Results for Pitch Mode (Test AFXZL1) 
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Figure 29. Improved Results for Pitch Mode (Test AFXZL2) 
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